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In this essay, we first sketch the development of ideas on the extraordinary 
dynamics of integrable classical nonlinear and quantum many body Hamil- 
tonians. In particular, we comment on the state of mathematical techniques 
available for analyzing their thermodynamic and dynamic properties. 

Then, we discuss the unconventional finite temperature transport of inte- 
grable systems using as example the classical Toda chain and the toy model 
of a quantum particle interacting with a fermionic bath in one dimension; we 
focus on the singular long time asymptotic of current- current correlations, we 
introduce the notion of the Drude weight and we emphasize the role played by 
conservation laws in establishing the ballistic character of transport in these 
systems. 
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1. INTRODUCTION 

The extraordinary stability of solitons upon collisions in integrable non- 
linear systems was first discovered in a numerical simulation of the Korteweg- 
de-Vries evolution equation, a system commonly studied in hydrodynamics 
and plasma physicsB. This discovery was soon followed by the development 
of a beautiful mathematical theory - the Inverse Scattering Method (ISM) - 
that allows the analytical evaluation of the time evolution of an initial pulse 
configuration using linear operations!; this is the analogue of the Fourier 
Transform for linear systems. These seminal works opened the new and still 
rapidly expanding field of nonlinear physics, with developments ranging from 
mathematical physics to applications^. 

Here, we should point out that there is a fundamental distinction be- 
tween integrable systems (mostly one dimensional) characterized by the pres- 
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ence of "mathematical solitons" and those with "topological solitons". On 
the one hand, the stability of mathematical solitons is guaranteed by a subtle 
interplay between dispersion and nonlinearity; this interplay is expressed by 
the existence of a macroscopic number of conservation laws constraining the 
dynamical evolution. On the other hand, the stability of topological solitons 
is enforced by a topological constraint; of course there are examples, like 
the sine-Gordon field theory, which possess topological solitons but they are 
also integrable. In this work we will focus on the transport properties of 
integrable systems. 

Due to the presence and stability upon scattering of nonlinear excita- 
tions, integrable systems are expected to show unconventional finite temper- 
ature transport properties, as ideal thermal, charge or spin conductivities, 
ballistic rather than diffusive transport. Within the traditional framework 
of linear response theory, the finite temperature dynamic correlations char- 
acterize the transport behavior and they are directly linked to experimental 
observations. However, although integrable models are considered as exactly 
soluble, meaning that the initial value problem can be exactly analyzed using 
the ISM, rather little progress has been done in the evaluation of dynamic 
correlations which remains at best technically very involved. 

In the quantum domain, parallel to developments on the analysis of 
classical integrable nonlinear systems, in the early 80's it was realized that 
the exact solution of a certain class of one dimensional quantum models 
by the Bethe ansatz (BA) method was equivalent to a quantum version of 
the ISM0. In this class belong well known prototype systems as the Hub- 
bard or spin-1/2 Heisenberg model, commonly used for the description of 
(quasi) one-dimensional electronic or magnetic materials. The BA method 
provides the exact eigenfunctions and eigenvalues and by a certain proce- 
dure (and assumptions) the exact thermodynamic properties and excitation 
spectrum. Similarly to their classical counterparts, the quantum systems 
possess a macroscopic number of conservation laws, characteristic of their 
integrability. It should therefore come as no surprise the proposition that 
quantum integrable systems should also exhibit unconventional transports. 
The situation however is similar to the one of classical systems; although 
the exact eigenfunctions and eigenvalues are known, the calculation of finite 
temperature dynamic correlations is still out of reach for most of the models. 

In the following, we discuss two simple examples, one classical and one 
quantum, in order to show the ideal conducting properties of integrable 
systems. Our strategy is to focus on the long time asymptotic value of the 
current correlations in order to demonstrate the ballistic transport instead 
of attempting to evaluate the full frequency dependence of the conductivity 
(or mobility). Furthermore, instead of linking the ideal conductivity to the 
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dynamics of soliton excitations, not very transparent for quantum many 
body systems, we directly relate it to the conservation laws characterizing 
integrable systems. 

First, we present a study of the energy current - current correlations 
(related to the thermal conductivity) for the classical Toda chain. Second, we 
introduce the notion of the Drude weight as a criterion of ideal conductivity 
and evaluate it exactly using the Bethe ansatz method in the context of the 
mobility for a toy model describing a quantum particle interacting with a 
fermionic bath. 



2. A classical system: the Toda chain 

The classical Todalattice is a prototype model for studying the physics 
of nonlinear excitations^. It is one of the first models analyzed using the In- 
verse Scattering methodH, the conservation laws characterizing this system! 
have been presented and it has even been invoked in attempts to describe 
nonlinear transport in DNA molecules!. As we mentioned above, although 
the initial value problem and the thermodynamic properties can be analyt- 
ically studied, there is no clear picture on the finite temperature dynamic 
correlations!^. 

A physical quantity of interest in an anharmonic chain is the heat con- 
ductivity. In a generic case, it is expected that the energy current correlations 
decay to zero at long times and the decay is fast enough so that a transport 
coefficient can be defined. This behavior seems rather difficult to observe in 
several one dimensional systems! 11 }! 12 !, where the currents decay to zero but 
often too slowly, leaving the issue of diffusive transport controversial. 

For an integrable model, ideal conducting behavior is expected with 
current correlations decaying to a finite value at long times. To quantify 
the contribution of nonlinear excitations, different ingenious methods have 
been devised! (soliton counting procedures). Here, we will use the long time 
asymptotic value of current correlations as a measure of ideal transport, 
related in integrable systems to the existence of conservation laws. 

To establish this relation, we will use an inequality proposed by MazuiO, 
linking the long time asymptotic of dynamic correlations functions to the 
presence of conservation laws: 

lira U\A(t)A)dt>j:^f- (1) 

T^oo 1 Jo - {Q^} 

Here () denotes thermodynamic average, the sum is over a set of conserved 
quantities Q n , orthogonal to each other {Q n Q m ) = (Qn)^n,m and we suppose 
that (A) = 0. 
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The classical Toda Hamiltonian for a chain of L sites with periodic 
boundary conditions is given in reduced units by: 

H = J2^+e- qi (2) 
1=1 z 

where pi is the momentum of particle I, xi its position and qi = xj+i — x±. 
The energy current for a system of interacting particles is given byllil: 

j =^ Pl h t + qie qi (3) 

1=1 1 

2 

where h t = \ + ^(e -9 ' + e~ qi - 1 ). 

We consider dynamic correlation functions in the fixed temperature- 
pressure thermodynamic ensemble: 

(A(t)A) = Z- 1 H d Pl d qi A(t)Ae-^ H+p V (4) 

2=1 

where Z = J ]\f = idpidqie~^ <yH+PL \ L = Ya=iQ.U P ls the pressure and f3 
the inverse of the temperature. 

In this thermodynamic ensemble, equal time correlation functions can 
be calculated analytically. For instance the average distance is given by: 

(q) = HP) - 9(J3P) (5) 

where ^f(z) is the digamma function. 

The classical Toda lattice is characterized by a macroscopic number of 
conservation laws. The first few ones are: 
L 

Qi = Em ( 6 ) 
i=i 

Qi = zbi + e- 91 (7) 

1=1 z 

Qz = zZ4 + (Pi+Pl+l) e ~ qi ( 8 ) 

1=1 6 

Qi = E^ + (^+m + i+pli)e-Hje- 2 He- ,! 'e-^ (9) 
i=i 4 1 

L 5 

= Y,^ + (Pi+PiPi+i+PiPf + i+Pi + i)e- qt (10) 

1=1 D 

+ (p, + pi + i)e- 2qi + (p, + 2 Pl+1 + W+2 )e-' ?i e- 9; + 1 ... (11) 
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with Q\ the total momentum, of course present in all translationally invari- 
ant systems, integrable or not, as also Q2 the total energy. According to the 
standard Green-Kubo formulation of transport theoryH-3 "subtracted fluxes" 
should be used in the dynamic correlation functions determining the trans- 
port coefficients. So in the case of energy transport we will study the decay 
of the "subtracted" energy currenttJ: 

~-E _ -E (Ql3 E ) n n0 ) 

j -j --^a (12) 

We see that the use of a subtracted flux is equivalent to removing the con- 
tribution of Qi in the long time asymptotic boun difor (j E (t)j E ). 

We will now calculate a bound on the long time asymptotic value of 
{j E {t)j E } by the Mazur inequality eq.([l]) using the first m conservation 
laws. We should note that Q3 has a structure very similar to the energy 
current, so we expect a large contribution from this term; actually in some 
quantum models like the spin-1/2 Heisenberg chain or the t-J model, the 
energy current is identical to a conservation law, directly implying a nonde- 
caying energy current and thus infinite thermal conductivity^. Here, Q' n s 
with n =even do not couple to j E so we will consider only Q n ,n = 3,5,7. 
Higher Q' n s can of course be included but the calculations become rather 
cumbersome. Orthogonalizing the conserved quantities which appear in the 
right hand side of eq.(|]) is equivalent to evaluating the expression: 

C™ E]E = Q E \Q)(Q\Q)- 1 {Q\j E ) (13) 

where (Q\Q) is the m x m overlap matrix of Q' n n and (Q\j E ) the overlap 
vector of j E with the Q' n s. 

~ 1 

In Fig.l we show the temperature dependence of C™^j E / (j E ) for m = 
1 (n = 3), m = 2 (rt = 3,5) and m = 3 (n = 3,5,7). At low T the 
behavior is linear with slope ^ for m = 2 and ^ for m = 3. It is interest- 
ing that this value is comparable to the value for the density of solitons0 
N s /N = ln(2)/-K 2 T. So, we find that the long time asymptotic value of the 
subtracted energy current correlations is finite and most interestingly that 
it increases with temperature. This trend we can interpret as evidence for 
an increasing contribution of thermally excited nonlinear excitations on the 
ballistic transport. 

The idea presented here provides, on the one hand a conceptual under- 
standing of the role played by the conservation laws on the finite temperature 
dynamic correlations and on the other hand a simple analytical method for 
evaluating, or at least giving bounds on their long time asymptotic value. 
A similar analysis can be carried out for quantum integrable systemsll3 al- 
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Fig. 1. Lower bounds on the long time asymptotic value of energy current 
correlations. 



though the complexity of the quantum conservation laws renders their wide 
use rather limited. 



3. A quantum system: the "heavy particle" model 

The Drude weight D (or charge stifness) was introduced as a criterion 
of an ideal conducting or insulating state at zero temperatureEl and recently 
extended as a measure of ideal transport at finite temperatures!. Within 
linear response theory, it is essentially the prefactor of the low frequency 
reactive part of the conductivity, a" = 2D a finite D implying a 
freely accelerating system. For a normal, diffusive system the Drude weight 
is zero at any finite temperature in the thermodynamic limit; according to 
the standard scenario the weight at zero frequency spreads to a "Drude 
peak" with width proportional to the inverse of the collision time. As we 
will see below, in integrable systems the Drude weight remains finite at all 
temperatures indicative of ballistic rather than diffusive transport. 

The Drude weight can be conveniently evaluated^ as the thermal av- 
erage of curvatures of energy levels e n of the system subject to a fictitious 
flux (j), 
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"-k^w^ (14) 

where p n are the Boltzmann weights and the sum is over all eigenstates 
of the system. It is also e^ual to the long time asymptotic value of the 
current-current correlation: 



if 



D = ^<i(t)i)k-*oo = ^£^IR,»| 2 (15) 

n 

so useful bounds on D can be obtained using the Mazur inequality following 
the same formulation as in the previous example. 

For integrable quantum many body models it can evaluated exactly fol- 
lowing recent developments in the Bethe ansatz technique, thus providing 
essential information on the transport properties of these systems without 
requiring the full calculation of the frequency dependence of the conduc- 
tivity. This type of analysis, still under discussion as it is technically in- 
volved, has been carried out for several one dimensional integrable quantum 
models as the Hubbard chain, the spin 1/2 Heisenberg and the nonlinear-u 
modelEiEiEi These calculations show that in most of the cases the Drude 
weight is finite at all temperatures implying ideal thermal, charge or spin con- 
ductivity. Recently, investigation of the finite temperature of these systems 
was also carried ojrtby a semiclassical approach^! and within the Luttinger 
liquid descriptioncJ'Ej. 

Here, we will demonstrate the idea behind this type of Bethe ansatz 
analysis by evaluating the Drude weight related to the mobility of a quantum 
particle interacting with a bath of fermions in a one dimensional svstemEa 
A similar analysis was carried out for a particle moving on a lattice^ but the 
case discussed below is simpler and shows a qualitatively new behavior. 

We consider a particle with coordinate y moving on a system of length L 
with periodic boundary conditions and interacting with a set of N fermions 
described by the coordinates 5— function interaction of strength c, 

^=-E^-S + 2cE^-y)- (16) 

j j 9 3 

When the mass of the "heavy particle" is equal to the mass of the 
fermions the model is integrable and so we expect a ballistic behavior of 
the mobility and therefore a finite Drude weight. To evaluate D using eq. 



(14), we consider the dependence of the energy levels on a flux eft acting 
only on the heavy particle. The momenta kj and the collective coordinate 
A describing the Bethe ansatz wavefunctions are then given by the following 
standard equations obtained by applying periodic boundary conditions, 
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Lkj = 2-kIj + 8(kj - A), j = 1, N + l, (17) 
0(p) = -2tan _1 (2p/c), (18) 

AT+l 7V+1 

Lj2 k i = 27T J2 1 i + 27rJ + ( 19 ) 

The eigenstates are characterized by the quantum numbers (Ij, J) and their 
energy is given by: 

N+l N+l 
E=Y.< k 3)=H k l (2°) 

3=1 3=1 

These equations can be solved to order 1/L as we consider the effect of the 
one particle on the ensemble of fermions. 

fc, = fcJ + ^(fc i -A), fc° = ^ (21) 

Thus the total energy can be written as, 

E = J2<k<j) + ^k°9(k°-A), (22) 

3 

l£*$-A) = ^ + *. (23) 
j 

Going to the continuum limit we obtain: 

E(p(k),A) = ±-J dkp{k){k 2 + jk6(k - A)), (24) 

P + <j)= — [ dkp(k)9(k- A), P = ^~. (25) 

Now we can define a correlation energy e c (A) assuming that the distri- 
bution of the fermion momenta is not affected by the presence of the extra 
particle and replacing the density p(k) by the Fermi-Dirac distribution f(k), 

e c (A) = — i dkf(k)2k8(k - A). (26) 
2ir J 



Using this formulation and the definition of the Drude weight eq.(|14|) we 
obtain: 



0.2- 



0.3- 



— c=2 

— c=4 

— c=8 



0.1 



0, 







1 



2 

T 



3 



4 



Fig. 2. Drude weight of the "heavy particle". 



where 



<9#(£; - A) 



(28) 




dAg(A)w(A), w(A) = e 



■/MA) 



(29) 



In Fig. 2 we show the normalized Drude weight D/Dq for different 
values of the interaction c as a function of temperature, with Dq = D(T = 



0) = (Tr/2)(tan- 1 (2k F /c)-{2k F /c)/{l+(2k F /c) 2 ))/tan- 1 (2k F /c) 2 and k F = 



Trn. The chemical potential is chosen so that we consider density n = 1; 
upon scaling n — > nc, (3 -> (3/c 2 the Drude weight D remains the same. 
We note that the behavior of D is not monotonic, initially decreasing at 
low temperatures because of the interaction and then tending to the free 
particle value at high temperatures. This is in contrast to the Drude weight 
of systems on a lattice (tight binding models) where D goes always to zero 
as (3 at high temperatures. The difference in behavior can be attributed to 
the bounded spectrum in lattice models in contrast to the unbounded one 
for continuous models. Furthermore, numerical calculations on this model 
show that, the Drude weight vanishes at any finite temperature when the 
mass of the "heavy particle" is not equal to that of fermions, as expected for 
any normal system. In conclusion, the presented analysis demonstrates the 
basic features of the generic finite temperature ballistic transport behavior 
of integrable quantum many body systems. 
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